Radiation from cosmic string standing waves 
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We have simulated large-amplitude standing waves on an Abelian-Higgs cosmic string in classical 
lattice field theory. The radiation rate falls exponentially with wavelength, as one would expect 
from the field profile around a gauge string. Our results agree with those of Shellard and Moore, 
but not those of Vincent, Antunes, and Hindmarsh. The radiation rate falls too rapidly to sustain 
a scaling solution via direct radiation of particles from string length. There is thus reason to doubt 
claims of strong constraints on cosmic string theories from cosmic ray observations. 
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Introduction Cosmic strings are one-dimensional 
topological defects which may have been created by a 
phase transition in the early universe |ll|. (For reviews 
see [§,|[ ■) As the universe evolves, intercommutations be- 
tween long strings produce oscillating loops. In the stan- 
dard scenario, these loops lose energy by gravitational 
radiation and eventually disappear. This produces a scal- 
ing solution where the average distance between strings 
is a constant fraction of the Hubble length. Most of the 
energy in the string network is emitted as gravitational 
waves, which we cannot observe, and only a small frac- 
tion appears as high-energy particles. 

However, in a recent paper B, Vincent, Antunes, and 
Hindmarsh claim that energy in a string network is lost 
by direct particle emission from long strings, rather than 
in gravitational waves. To back up this claim they study 
large-amplitude sinusoidal standing waves, and claim 
that the energy emission rate is sufficient to explain 
scaling behavior with the great majority of the energy 
emitted as particles. Moore and Shellard S found that 
the emission rate fell exponentially with wavelength, but 
their amplitudes g were much less than those of [Q. 
Furthermore, the range of wavelengths in which [Q saw 
exponential fall off had no overlap with the wavelengths 
studied in ||] . 

Here we simulate the same large-amplitude waves as in 
yi and cover part of the same wavelength range, but we 
come to a different conclusion. In our simulations, the 
energy emission rate declines exponentially with wave- 
length, and thus cannot account for the large direct- 
emission rate claimed in 0] . 

Model We work with the Abelian-Higgs model, which 
produces local strings with no massless degrees of free- 
dom in the vacuum. The Lagrangian is 




FIG. 1. A surface of constant energy density shows a sinu- 
soidal wave on a cosmic string and the energy emitted from 
the high curvature regions. 
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We work with units such that rj = 1 and e = 1, and we use 
the "critical coupling" regime in which j3 — A/(2e^) — 1 
so that in our units A = 2. As in Q, we study strings 
whose initial core position is given by a sinusoidal wave 
y — Acoskx, with amplitude A — A/2 = n/k |0]. 

A preliminary investigation shows that emission from a 
standing wave is not uniform in time, but rather consists 
of a series of bursts emitted when the string is momentar- 
ily stationary with large amplitude waves in its position. 
Figure |] shows a snapshot of the energy density around 
a string at one such point, and Fig. || shows a plot of the 
energy emitted over several oscillations. 

It thus appears that standing wave radiation is akin to 
radiation from cusps, and results from the overlap of the 
tails of the the string fields. We use such a model below 
to compute a theoretical expectation of the dependence 
of radiation rate on wavelength. 

Expectations Vincent, Antunes, and Hindmarsh B 
argued as follows: in a scaling network, the distance be- 
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FIG. 2. Energy emitted from one wavelength of a standing 
wave. The radiation appears as a series of bursts. 



tween strings, £,, scales with the Hubble distance, which is 
proportional to time. In a volume ^'^ there will be string 
length roughly ^, so the energy density in the string net- 
work is p = IJ./S,'^, and thus p = —2fj,^/S,^. Since ^ is a 
constant, p^^ is constant. As a model they used a si- 
nusoidal standing wave with wavelength A in a box of 
volume A^. They expected pX^ to be constant. If we let 
E be the energy of a single wavelength of the string, then 
E = pA^, and thus E is independent of A. If we let Pl 
be the power per unit length radiated from the standing 
wave, we need 
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to sustain a scaling network from energy emission of this 
type. 

In contrast, analyzing the fields around the string 
would lead to a different conclusion. A straight, static 
string is topologically stabilized in a minimum-energy 
configuration, and so cannot radiate. If the string is 
curved, then there is the possibility for radiation, but 
since the fields fall off exponentially toward the vacuum 
at large distances from the string, one would expect the 
amount of radiation to be suppressed by an exponential 
factor depending on the radius of curvature, R. This 
seems in keeping with Fig. H^, which shows the radiation 
coming from the points of maximum curvature. 

As a specific model, one can imagine that an element of 
momentarily stationary curved string gives up an amount 
of energy proportional to exp{—aR)dl, where a is a con- 
stant of 0(1) and dl is the length of the element of string. 
The total energy emission is then 
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For a sinusoidal wave, y = A cos kx, the radius of cur- 
vature is 
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We will consider the region around one of the peaks of 
the sinusoid. The energy emission is dominated by the 
region where x is near zero, so that R is small. In this 
regime, we can approximate 

jl+A^k^X^)^/^ 1 /3 2 1\ 2 .^ 

In our case, we are going to consider a fixed ratio of 
amplitude to wavelength, A = A/2 ~ n/k so we get 
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We can now do the integral of Eq. (W, approximating 
dl = vl + Ti'sin kx dx by just dx, and extending the 
limits of integration to infinity, to get 
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Since we keep the amplitude a fixed multiple of the 
wavelength, the period of the standing wave is just pro- 
portional to A. If we consider a half wavelength of string, 
it emits bursts of energy E twice per cycle, so the power 
is ex X~^''^e~^^, and the power per unit length is 
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Simulation The simulation is based on a lattice ac- 
tion, as described in Iq]. However, in the present case 
we have used a different lattice spacings in the 3 cardi- 
nal directions. The maximum speed of the string is quite 
large, and leads to a Lorentz contraction of the field pro- 
file in the direction of motion. To accurately represent 
the fields, the lattice spacing should be proportional to 
l/7i — -^/l — vf, where Vi is the component of the string 
velocity along axis i. In the z direction, where there is 
no motion, we have used a lattice spacing of 0.33, which 
seems to be the largest that gives reliable results. The 
corresponding spacings in the x and y directions are 0.31 
and 0.10 respectively. The Courant condition requires 
At < (Ax"" + Aj/~" + Az"")~^^^ « 0.09, and in our 
simulations we use At = 0.08. 

To extract the energy which is emitted by the string, 
we have used absorbing boundary conditions on the y 
and z faces and accumulated at each step the amount of 
energy that they absorb. The conditions are 



Ey = n X B . 
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where D denotes the covariant derivative, n the out- 
ward normal unit vector on each boundary, and E^ = 
E — n(E • n) the transverse component of the electric 
field. This corresponds to the zeroth-order absorbing 



boundary condition for free electromagnetism g. The 
energy flowing into the boundary is 

Eabsorbed — / S • dn (10) 

JQ. 

where 57 is the boundary surface, and S is the Poynting 
vector, given by 

Sj = -Dt^DjCJ) - Dt(f>Dj4> + (E X B)j . (11) 

Using Eqs. (H), we can rewrite Eq. (nffl) as 
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which is easily computed. 

To produce the sinusoidal waves, we use the same tech- 
nique that we used to generate cusps in g, i.e., we cre- 
ate two traveling wiggles on a straight string that will 
combine to produce the desired sinusoidal form. (In 
the Nambu-Goto approximation, the resulting sinusoid 
would be exact; in our case there will be a distortion 
of the shape due to the dynamics that occur before it 
is formed, but this effect will be small because the wig- 
gles out of which the sinusoid forms are not themselves 
strongly curved.) The initial field configuration for a 
moving wiggle is known exactly, from a result of Vachas- 
pati [|l^ . The advantage of this technique is that it does 
not require the use of relaxation, as is necessary in other 
field theory simulation schemes MM- 

The two original wiggles will overlap to form a single 
wavelength of the standing wave, from one minimum of 
y to the next. At this point it is possible to change to 
periodic boundary conditions in the x direction, so that 
the straight part of the string is removed, and we are 
left with a single wavelength of standing wave in a pe- 
riodic box. This technique was used to produce Figs. |l| 
and y, but it is not an accurate method for extracting 
the radiation rate, because energy coming from differ- 
ent burst is not clearly separated by the time it reaches 
the boundaries. When each burst of radiation is emit- 
ted, the string changes its shape and amplitude and its 
subsequent evolution does not correspond to a constant 
amplitude sinusoid any more. Of course, for large enough 
A this would not matter, but for the wavelengths in the 
range of our simulation, it makes a significant difference. 

To avoid this problem, we allow the original wiggles 
to pass by each other beyond the point of the overlap, 
so that they generate just a single burst, and then sep- 
arate. The place from which the burst is emitted is at 
the center of the overlapping region, and the string near 
that point has been following the same evolution as in a 
real standing wave for a half period, so we feel that this 
burst accurately represents a single burst of a standing 
wave oscillation. 

Using the expressions given above for the Poynting 
vector on the boundaries, we can compute the energy 




FIG. 3. Absorption of a burst of energy. The two wig- 
gles have passed through each other and are now separat- 
ing. Three sides of the bounding box are shown, with darker 
shades showing areas of greater energy absorption. 



absorbed at each time step in our simulation. The en- 
ergy absorption on the box faces is shown in Fig. ^. 
Integrating all the energy absorbed from the first burst 
of radiation we can compute the power emitted by a si- 
nusoidal standing wave. 

We have repeated this procedure for different values 
of the wavelength ranging from A = 12 to 30, in natural 
units. In Fig. ^ we plot the power emitted per unit length 
and compare it with theoretical predictions from Eqs. (|^) 
and (H). We see that the exponentially decaying model 
fits quite well. The line shown has the form 
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with _Ro ^ A/(27r^), but we do not have sufficient accu- 
racy to confirm the exponent of A or the exact constant. 
A curve with X~^ or A"^ and somewhat different constant 
would fit equally well. On the other hand, the form of 
Q does not fit at all. 

Discussion We have simulated large-amplitude 
standing waves on local cosmic strings, and found an 
exponential decrease in radiation power with increasing 
wavelength. Our technique proceeds from separated wig- 
gles with exact initial conditions, and we have used quite 
a small lattice spacing as compared to other authors, so 
we feel that our results are reliable. 

The constant in the exponential gets its dimensions 
from the string thickness (10~'^''cm for a GUT scale 
string), and thus for waves of any reasonable cosmologi- 
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FIG. 4. Power per unit length from simulation (points) 
compared with exponentially decaying model (solid line) and 
linearly decaying model of li (dotted line). 



cal size, the radiation is utterly negligible. One could in 
principle imagine that strings in cosmological networks 
still have excitations at wavelengths comparable to their 
thickness, but this does not seem reasonable. Such wig- 
gles will be rapidly smoothed out by gravitational radia- 
tion, and there is no mechanism for regenerating them at 
such small scales. Thus we conclude that direct radiation 
of particles from string length cannot play a significant 
role in the production of cosmic rays or the maintenance 
of a scaling network. As a result, cosmic ray observations 
do not rule out field theories that admit cosmic strings, 
as claimed in |,|ll|,|l|. 
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